Dynamic stability of crack fronts: Out-of-plane corrugations 
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The dynamics and stability of brittle cracks are not yet fully understood. Here we use the Willis- 
Movchan 3D linear perturbation formalism [J. Mech. Phys. Solids 45, 591 (1997)] to study the 
out-of-plane stability of planar crack fronts in the framework of linear elastic fracture mechanics. 
We discuss a minimal scenario in which linearly unstable crack front corrugations might emerge 
above a critical front propagation speed. We calculate this speed as a function of Poisson's ratio 
and show that corrugations propagate along the crack front at nearly the Rayleigh wave-speed. 
Finally, we hypothesize about a possible relation between such corrugations and the long-standing 
problem of crack branching. 



Crack propagation is a major vehicle for material fail- 
ure 0, Q. Yet, a complete theoretical understanding of 
the dynamics and stability of rapidly propagating tensile 
cracks in brittle materials is lacking [!, H[. In the last 
few decades it has been experimentally established that 
rapid cracks undergo various types of dynamic instabil- 
ities. Planar crack fronts can experience a 3D micro- 
branching instability [sHl], a 2D oscillatory instability 
distributed micro-cracking in their vicinity [UO 
and support propagating 3D waves, with both in-plane 
and out-of-plane components 14H17|. 

Recently, progress in understanding these phenomena 
has been made. For example, linear in- plan e crack front 
waves have been theoretically predicted |18M2l| , dynamic 
crack branc hing has been observed in 2D numerical sim- 
ulations U 25j . en erg etic bounds on 2D crack branching 
have been derived 26l430| and the 2D oscillatory insta- 
bility has been explained based on intrinsic nonlinear- 
ities (3ll [lH , whose role in dynamic fracture has been 
highlighted in 33-3^]. While progress in understanding 
3D crack front instabilities in the context of quasi-static 
mixed- mode fracture 3, 41 1 and dynamic tensile frac- 
ture [42| has been made very recently, the 3D out-of-plane 
nature of crack front waves [3| and the micro-branching 
instability Q remain theoretically elusive. 

In this Letter we study the out-of-plane stability of pla- 
nar crack fronts in the framework of linear elastic fracture 
mechanics using the Willis-Movchan 3D linear perturba- 
tion formalism The emergence of out-of-plane 
crack corrugations within a minimal scenario is discussed. 
We calculate the critical propagation speed for the onset 
of such corrugations, which propagate along the crack 
front at nearly the Rayleigh wave-speed, as a function 
of Poisson's ratio. We then speculate about a possible 
relation between these corrugations and the 3D micro- 
branching instability and its fractographic consequences. 
Our work is directly inspired by [42| . 

To set the stage for our stability analysis, consider a 



planar crack propagating steadily in an isotropic linear 
elastic solid and introduce small perturbations which give 
rise to the following perturbed crack surface configura- 
tion 

S e (t) = {x:~co<xi<vt + e<j)(x2,t), 
-oo<:i'2<oo, x 3 = eyb(xi,X2,t)} , (1) 

where {x\, xi, £3} is a fixed Cartesian coordinate system 
and v is a constant propagation speed in the xi-direction, 
smaller than the Rayleigh wave-speed cr. The smooth 
and bounded functions 4> and ip describe in-plane and 
out-of-plane perturbations, respectively, and < e <C 1 is 
a small dimcnsionlcss amplitude. Note that in Eq. ([1]) 
we use 00 to schematically represent finite macroscopic 
length scales that are much larger than any length scale 
characterizing the perturbations. Such scales will be ex- 
plicitly invoked below and will play a role in our analysis. 

The isotropic linear elastic solid is characterized by a 
stress field cr and a displacement field u that satisfy the 
momentum balance equation V • cr = pdttu, where p is 
the mass density, and are related through Hooke's law 
cr — Atre J + 2/ie, where e= ^[Vit + (Vw) T ]. Here J is 
the identity tensor, and A and p are the Lame constants 
that are related to the dilatational and shear wave-speeds 
\/(A + 2/i) / p and c s = yj ' p/ ' p, respectively. 

0, where i = 

1—3 and n is an outward unit normal to S e (t), and some 
external conditions are specified on the outer boundaries 
of the solid. This formulation constitutes Linear Elastic 
Fracture Mechanics (LEFM) Q. 

A fundamental prediction of this framework is that the 
asymptotic crack front expansion of the stress field cr is 
dominated by a universal square-root singularity in the 
distance from the front. For the unperturbed crack (Eq. 
(fTj) with e = 0), under pure (symmetric) tensile loading 
conditions, the expansion of the tensile stress component 
ahead of the front takes the form [l[ 
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where X = x± — vt. Here K^ is the mode I (tensile) 
Stress Intensity Factor (SIF), a fundamental quantity 
that quantifies the intensity of the linear elastic singu- 
larity, is the coefficient of the sub-leading term to 
be used below and the superscript (0) refers to the un- 
perturbed crack. Note that Kj ^ and are assumed 
to be independent of x 2 and t. The crack front perturba- 
tions described in Eq. ((T|) break both the translational 
symmetry along X2 and the x 3 — > — x 3 symmetry of the 
global loading. Therefore, the stress er associated with 
the perturbed crack admits the same asymptotic expan- 
sion (ahead of the front, along its local tangent) as in Eq. 
([2|), but now Ki{x2,t) depends on X2 and t, and there 
exist also mode II (in-plane shear) and mode III (tear- 
ing) singular contributions proportional to Ku{x2, t) and 
Kiu(x2,t), respectively. 

The major goal of a linear perturbation theory is to cal- 
culate the SIF's Ki 11 in to linear order in e: Kj(x 2 ,t) = 



that K^-a^VL and A 3 V > = Kf ! a 3 (v)/L, where a 3 (v) 
is a real dimcnsionless function. Furthermore, the Fourier 
transform of Tn(x2, t) can be expressed as 

Tu (6 , w) = K ( p In (w/ LW(& , w) , (7) 

where tu is a complex function. 

Therefore, the Fourier space version of Eqs. ©-(HI 
reads 

#$(6,00= [-|6|ei33ii(w/|6l) + *(w/«)wi3 (8) 

-Giaos/i + VM/Lhdoj/ |6| , 161 L)]r(&oj)Kf } , 

0(6 ,o;)#f, (9) 
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K^ ) 6 j i + eK^\x 2 ,t), with j = 1,11,111. This problem 
was solved by Willis and Movchan in 1997, yielding [HI 
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K { I 1 \x2,t)=Q 33 ^K i ° ) + v^^A 3 0) , 
*r#(a; 2 ,t) = [-GiaQn* ^-wiaW]^ (4) 
- O 13 ii*A 3 0) +T n (x 2 ,t) , 

where * stands for a convolution in (x 2 , £) and ?/>* is ip 
evaluated at the crack front. We do not present here 
Km because it plays no role in what follows. 

Let us discuss the physical meaning of the various 
terms in Eqs. ([3|)-(j4]). Qij are weight-functions that 
quantify the non-local (both in space and time) elasto- 
dynamic interaction of the crack front with itself. They 
admit explicit Fourier transform representations Q,-., and 
are homogeneous functions of degree one (see [42|, EH - li^ 
for details). Tu(x2,t) involves convolutions of additional 
weight-functions with the "history" of out-of-plane per- 
turbations, i.e. xj) for X < 0, and effective tractions asso- 
ciated with zero order stresses [ijj . The other terms 
are purely geometric in nature and correspond to local 
rotations and translations of the crack front. Finally, the 
kinematic functions O13 (v) and ui 3 (v) take the form 



13 (v) = -2a{/3 - a)(l + P 2 )/R{v) , 



(5) 



wis(w) = -[8aP - (1 + P 2 )(2 + a 2 + af3)]/R(v) , (6) 

where a 2 = 1 - v 2 j<? d , (3 2 = 1 - v 2 /c 2 and R(v) =4a(3 - 
(1 + (3 2 ) 2 (whose root is the Rayleigh wave-speed cr). 

In order to prepare Eqs. ©-(jl]) for the stability analy- 
sis to follow, we first Fourier transform all relevant func- 
tions with respect to (x2,t) such that /(6i w ) is the 
transform of f(x2,t). Noting that to leading order in e 
one has d Xl tp* = w _1 9t'0*, the wave-number of the pertur- 
bation ip in the propagation direction xi is readily given 
by 6 = w/u. In a minimal scenario, a crack problem 
features one macroscopic length scale L and one macro- 
scopic stress scale a°° . Dimensional analysis then implies 
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where we used the homogeneity property of Q to de- 
fine Q lli33 (6, u) = |6I 811,33 (w/ |6|). Equations 
which serve as a basis for our stability analysis, show that 
in-plane and out-of-plane perturbations are decoupled to 
first order in e [45| , though experiments demonstrate such 



a coupling [1J] , suggesting the intervention of higher or- 
der corrections. 

In order to study the out-of-plane stability of the crack, 
we need an equation of motion for the front in terms of 
the SIF's. It is well-established that under quasi-static 
and 2D (i.e. when the front can be treated as a tip) con- 
ditions, the Principle of Local Symmetry (PLS) Kn = 



is valid Ml 51 



This implies that no intrinsic non- 
linear scales near the crack front play a role in crack 
instabilities and hence that crack patterns are always de- 
termined by extrinsic (typically geometric) length scales. 
However, it was recently shown that for strongly dynamic 
conditions (i.e. v of the order of cr) in 2D, a generaliza- 
tion of the PLS is required and intrinsic nonlinear scales 
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play a decisive role in rapid crack instabilities 
The role of intrinsic length scales in the context of quasi- 
static 3D crack propagation under mixed-mode Kj+Krrr 
conditions has also been demonstrated recently [40, |41( . 

Here, however, we seek a minimal scenario for the 
emergence of out-of-plane instabilities and hence wish 
to deviate as little as possible from LEFM. Therefore, 
we adopt the PLS, Ku = 0, even though we consider 
3D cracks under dynamic conditions. Finally, as a crack 
propagation condition we use the generalized Griffith cri- 
terion, quantifying the balance between energy flowing 
into the crack front region and the dissipation associated 
with the fracture process [l[, which together with the 



PLS reads [52| 



a{l - fi 2 )K 2 i/[2^R{v)] = T and K 



11 







(10) 



Here T is the fracture energy, assumed to be v- 
independent. Note that in general, the energy balance 
relation contains also terms proportional to K u and K\ n 
[H ; the former is omitted due to the Ku = condition 
and the latter because it docs not contribute to order e. 
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Before discussing the out-of-plane stability, we br iefly 
consider the in- plane dynamics of the crack front [19J. 
Expanding the energy balance relation in Eq. (|10[) and 

using K^\^ 2 ,tj) of Eq. ((9]), we obtain 



ail 



R(v) 



2ira,3 



. (11) 



In the limit | ^2 1 ^ 3> 1 this equation admits a real solu- 
tion for every v, giving rise to a dispersion relation of 
the form w/l&l = Sf(v) + 0(1/\&\L), where c 2 — sj + v 2 
is the propagation speed in a fixed coordinate system 
(sf is the propagation speed seen by an observer moving 
with the crack front). These are the well-known in-plane 
crack front waves whose speed c/ is plotted in 

Fig. QJ>. Note that while the linear perturbation theory 
determines a dispersion relation in terms of u and 161, 
it cannot select them separately. Moreover, although in- 
plane crack fronts might play a role in fracture dynam- 
ics, they cannot leave fractographic marks by themselves, 
which entails an out-of-plane component. 

To study the out-of-plane stability, we apply the PLS 
to Eq. ((SJ, i.e. set K^j =0. The outcome is a complex 
equation for a complex variable o->/|6l and a real variable 
\&\L. We expect Kjj =0 to be equivalent to 

9t[u/M=h(v,\^\L), 9?K|6|] = <?W6|£) , (12) 

where h(-) and g(-) are in principle calculable. The crack 
stability against out-of-plane perturbations is determined 
by the sign of 53[o;/|6|]- For 3[w/|6|] <0 perturbations 
are attenuated and the crack is stable. For 9[w/|6|] >0 
perturbations are amplified and the crack is unstable, de- 
veloping surface corrugations. Therefore, the instability 
threshold is determined by S*[cj/|6|] = EH- Be- 
fore we perform this analysis, let us discuss two possible 
scenarios. In the first, we set |£ 2 |£ — > 00. Therefore, 
at the instability threshold, the second equation in (fT2)) 
reads g(v,oo) = 0. The solution, which determines the 
critical speed v c for instability, is then substituted into 
the first equation in (fT2"j) , yielding !ft[w/|6|] = h(v c ,oo), 
which is a dispersion relation for corrugation waves (sim- 
ilar to the one derived for in-plane front waves). In the 
second scenario, we keep \&\L finite (albeit large) and 
obtain g(v, \&\L) = and 5R[w/|6] = h(v, \&\L). Here 
both ui and |6| are being selected, and the critical speed 
V c (if exists) is determined from a solvability condition. 
In this case, the wave-number in the propagation direc- 
tion, £1 =uj/v, is also selected and the complete crack 
surface topography is determined. 

To test these scenarios we use the stability condition 
3[w/|6|] = to obtain (for v<c R ) [H| 



23 [?ii(«)] = s d v ln[/3(l - (3 2 )/R(v)] 



(13) 



K)/=0 (see Eq. (©), we obtain 

w X3 = -v 3 [hi (s, |6l L)\ /(s VML) 



(14) 
(15) 



where s = w/|6l is now real. Substituting this result in 



where 2uj 13 = 2wi 3 - Q13V d v ln[/3(l - f3 2 )/R(v)]. 

To analyze these equations within the first scenario we 
set l&lL — > 00 such that Eq. (fl"5|) becomes 0)13(1;) = 0. 
This equation admits no solution as ojia(v) <0 for every 
v<cr. Indeed, it was recently demonstrated that within 
this scenario there exist no solutions with real w/161 an d 
that out-of-plane perturbations experience attenuation 
for every v [42[. To test the second scenario, we pro- 
pose an iterative scheme for solving Eqs. (|T4|) - (fl~5)) . As 
macroscopic length scales are assumed to be much larger 
than any length scale characterizing the perturbations, to 
leading order we can set l^al-^ - >oo on the right-hand-side 
of Eq. (fl"4"]). obtaining 3? [q n (s)] ~ 0. If the latter admits 
a solution, then s is determined to order 0(1/^/161 Li) 
and Eq. ([TSJ yields y/\&\L = -v$s [t H (s,oo)] /(su>ia), 
which determines |6I and requires the calculation of In. 
Once tn is known, the steps can be iterated to produce 
higher order corrections in powers of 1/ y|6l L. Here, 
however, we do not calculate tii, which is a complex 
function that depends on the geometry and loading con- 
ditions of any particular crack problem, and hence trun- 
cate the iterative scheme at its first step, 3?[g 11 (s)] ~ 0, 
assuming that the next steps yield a solution for |6|- 

To explore this scenario we look for solutions of 
^[<7ii( s )] = Oj which we denote by s c . It turns out 
that a real solution s c exists only above a critical prop- 
agation speed v c , which is a function of Poisson's ratio 
v= A/[2(A+/i)]. Examples are shown in Fig. QJi. Once s c 
is known for a given v, the crack corrugation wave-speed 
c c , in a fixed coordinate system, is given by c 2 . =s c + v 2 . 
In Fig. [TJd c c /cji is plotted vs. v/cr for two values of v 
and v>v c . c c is nearly the Raylcigh wave-speed cr. For 
completeness, we plot also the in-plane crack front wave- 
speed Cf, which exists for every v. In Fig. [5] we show the 
critical speed for the onset of out-of-plane corrugations 
v c as a function of Poisson's ration v. This constitutes 
our main result: the existence of an out-of-plane crack 
instability above a critical speed within a minimal sce- 
nario, accompanied by the propagation of corrugation 
waves along the crack front. v c is a monotonically de- 
creasing function of v, which vanishes above i/ c ~0.4. 

Might the out-of-plane corrugations discussed here be 
related to the micro-branching instability? Providing a 
definite answer goes well beyond the scope of the present 
work. It is, however, conceivable that the two phenomena 
are related. A micro-branching event, like the corruga- 
tions discussed above, breaks the translational symmetry 
along X2 (it is localized along this dimension (jj) and goes 
out-of-plane. One may then speculate that the out-of- 
plane front instability may trigger the micro-branching 
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FIG. 2. (Color online) The normalized critical speed v c /cr for 
the onset of out-of-plane instability (solid blue line) and the 
normalized critical speed Vb/cR above which 2D branching is 
energetically possible vs. v (dashed red line) [3(il ]. 



FIG. 1. (Color online) (a) 5ft [<j n ] vs. s/c a with ^ = 0.25 for 
v = Q.3cr (solid line, no zero-crossing) and v — 0.7cr (dashed 
line, zero-crossing), (b) The normalized crack corrugation 
wave-speed c c /cr vs. v/cr for v = 0.25 (solid blue line) and 
v = 0A (dashed blue line). For completeness, we show also the 
normalized in-plane crack front wave-speed c/ /cr for v — 0.25 
(dotted red line) and v — OA (dash-dotted red line). 



instability, or if the two instabilities are of different ori- 
gins, that one may amplify/suppress the other [53]. 

To see what the implications of this speculation might 
be, we need an estimate for the critical speed for micro- 
branching Vf,. In the absence of a better alternative, we 
use the energetic bound for 2D branching [3(|, which is 
shown in Fig. [2] together with v c - The two curves inter- 
sect at i/~0.32. Consider materials with v < 0.32 (e.g. 
i/~0.2 for silica glass) where Vf,<v c . For v<Vb<v c out- 
of-plane perturbations are stable and branching is not 
possible. For vi, < v < v c out-of-plane perturbations are 
still stable, though branching is energetically possible. 
Finally, for Vb<v c <v out-of-plane perturbations are un- 
stable and branching is energetically possible. Consider 
then materials with ^>0.32 (e.g. i/~0.35 for Plexiglas or 
v~ 0.5 for elastomer gels). Here v c < Vb and there exists a 
regime with v c < v < Vb, where out-of-plane perturbations 
are unstable, while branching is energetically prohibited. 

As corrugations and branching events leave traces on 
fracture surfaces, the different v regimes described above 
might be related to the common mirror-mist-hackle tran- 
sition, a fractographic characterization that corresponds 
to progressive increase in surface roughness with increas- 
ing crack speed [54 |. For example, in a regime where 
the crack front is stable against both corrugations and 
branching, we expect mirror-like surfaces. For v > v c ~ 
0.4, where such a regime does not exist, we do not expect 
mirror-like surfaces even at very low propagation speeds. 



In a regime where only one of the instabilities can be 
triggered, typically at higher u's, surface roughness will 
emerge, which might correspond to the mist zone. When 
both instabilities can take place, at yet higher u's, we 
expect denser surface features and enhanced roughness, 
which might correspond to the hackle zone. We hope to 
further explore these issues in the future. 

While our analysis, assuming LEFM and the PLS, im- 
plies that the scale of out-of-plane corrugations is deter- 
mined by an extrinsic (geometric or loading-inherited) 
length scale L, in practice such corrugations typically 
feature much smaller scales 0, Hij]. Nevertheless, we be- 
lieve our analysis remains relevant also when intrinsic 
length scales, associated with near crack front nonlin- 
caritics, microstructures and dissipation, play a role in 
selecting the scales of corrugations. The idea is that the 
SIF's uniquely couple the large scales of a crack problem 
to the small scales near the front. Therefore, whenever 
a region where the square-root singular SIF's fields dom- 
inate the mechanical response exists, which is the case 
for a broad range of materials the calculation of the 
SIF's in Eqs. (JSj)- is required. Moreover, we expect 
any generalization of the PLS - incorporating the role of 
intrinsic scales - to be represented in terms similar to tu 
and thus the structure of Eqs. (|T4")) - (|15|) to remain valid, 
where the role of L is played now by an intrinsic length 
scale. 

To conclude, we demonstrate - within a minimal lin- 
ear elastic fracture mechanics scenario - the existence of 
an out-of-plane crack front instability characterized by a 
Poisson's ratio dependent critical propagation speed. We 
speculate on a possible relation between this instability 
and the 3D micro-branching instability, and its fracto- 
graphic implications. Extending our anal ysis to explic- 
itly include intrinsic nonlinear effects j3ll l32| and finite 
geometric scales is an important direction for future in- 
vestigation. 
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Supplementary Information 

The aim of this document is to provide additional technical details about the calculations of the complex functions 
q 33 (s) in Eq. (11) and <7ii( s ) m Eqs. (13)-(14) in the manuscript. 

FIRST ORDER PERTURBATIONS OF THE STRESS INTENSITY FACTOR 

Willis and Movchan introduced a convolution identity (Eq. (3.12) of [Sl[) that allows to calculate the change of 
the stress intensity factors due to a change in the crack surface, induced by a perturbation of its propagation front. 
Starting from an unperturbed state defined by a large crack submitted to a mode I loading, the first order perturbation 
terms of the stress intensity factors, Kj (£2, w) and Kjj (£2,1*1) are given by (Eqs. (3)-(4) of the main text) 



Q 33 (£ 2 ,oj)K\ 



(0) 



ZL A (o) 



<M6,w) , 



K$(M = \-®i3Q n (&,u)Ki 0) +i(ui/v)ui 13 K^ 



e 13 4 0) 



</> (&,w)+Ti / (£ 2) w) 



(Ai6) 

(A17) 



Explicit expressions for the quantities in Eqs. (|A16[) - (|A17[) can be found in [S2J. The "history" term Tjj includes the 
contributions of the zero order traction terms on the broken surface convoluted with the corrugations of the crack 



surface ip(X < 0,X2,t). It is expressed as [S2 



T u (x 2 ,t) = [U] u * (PP) + [U} 21 * (P 2 (1) ) - (U) 31 * [P 3 (1) ] , 



(A18) 



where [•] and (•) denote the jump and the avera ge o f the corresponding quantities through the unperturbed crack 
plane. Here U are the weight functions defined in [S2J and P' 1 ' are traction terms obtained by imposing the boundary 
conditions on the edges of the perturbed crack. The latter are given by 



P. (1) = -d 



(A19) 



Also, one can show that Kj does not include such a term, i.e. that Tj(x 2 ,t) = 0. 
The functions -,(£ 2 , u>) are homogeneous of degree one and can be written as 



Qii(6,w) = |6l5iiH 161), 



?33 (M = |6 I 533 W 161) 



(A20) 



where [^2 lis the modulus of the complex wavenumber £2 and w is a complex frequency. The functions q xl and q 33 are 
given in [S2| 



111 S = "i s + ^ b (s) + -± -r±- 

% 3 (s) = -i&(s) + ka(*)+ I filn(T( S ,a,l)) 

I J c ZlT 



2tt 



ln(T(s,a,l)) , 



(A21) 
(A22) 



where 



T(s,£i,l) 



7 2 [4gj - <i) 2 /c 2 - gYg - gg/g - £1 + ((s - <!) 2 /^ - 2£\f 



COO 



6+ = 



T+(«,i,l) 



c2iZ(u)(-i -£+(«)) 



/3ir 



r_(a,-i,l) 



b- = M-i-tf(s) , 



(A23) 
(A24) 

(A25) 



with £ 2 = £ 2 + 1 and 



T-t(s, ±i, 1) = exp 



2ni 



ln(T( 8 ,fr,l) 
6 



(A26) 
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The contours C± encircle the branch points £^ and £^ respectively, with 



^ = ~T^f ±iqa] qa = lf--(£^' (A27) 



^ = -(^F ±%; 9c = W 1_ (^' (A29) 

with a 2 = 1 — v 2 /c 2 d , f3 2 = 1 — v 2 /c 2 , j 2 — 1 — v 2 jc\. The kinematic functions 0i3(?;) and wi3(i>) are defined in the 
main text. 

Finally, using the homogeneity of the functions U and P^ 1 ' and the time-independance of the zeroth order stresses, 
one can show that the function T/j can be written in Fourier space as 

T n fa,u)=K? ) ^hiM\b\Ab\L)r(.M , (A30) 

where tn is a complex function of degree zero and L is the geometrical length scale introduced in the main text. Here 
we used the fact that to leading order in e one has tp(X, X2,t) = ip(X + vt, x%). 



CALCULATIONS 



For a real s and V ' s 2 + v 2 < c R , Eqs. (|A21|) - (|A22j) reduce to 



where 



„ r , 1 ( a+a- — b+b- \ „ . , . , , . , 

M = -9c + x% + + + , + » F+ * , A31 
2 \a+a- + o + o_ / 

3 [^ S )]=(^lW +3[f+(s)1 ' (A32) 

» [933 (*)] = -Qc + \<la + » [f+ (*)] , ( A33) 

» [93s(*)] - T^W - ir 1 ^ + 9 [*+(*)] > (A34) 



F+(*)= [ ^HT(s,^l)), (A35) 



The functions a+a- and are given by 



6+6- = 1 + y/1 -H 2 /c 2 , (A36) 
R(v) (l + VI - #74 



c 2 i?(z;)(l + v/l-ff 2 /^) 
o+a- = a exp[G + (.s) + GL(*)] , (A37) 



with 

# 2 = s 2 + w 2 , (A38) 

and 
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The main calculations consist in computing the integrals F + (s) and G±(s). The algebraic manipulation of F + (s) was 
performed in S3|, |S4| and is given by 



F + (s) 



1 

2^ 



y/J{J - H 2 ){J -v 2 f 



(2 - J/c 2 ) 2 



dJ 



TT J C 2 (J-V 2 )' 



tan 



4^T—j/? d ^j/c[^T 

(2- J/c 2 ) 2 



(A40) 



It is now established that the energy balance of the in-plane crack front dynamics leads to non-dispersive waves with 
real velocity s c for all crack front propagation speeds [S4|. Using Eq. (11) of the main text, one then has 



®{q 33 (s c )}=Q 

[2q 33 (s c ) - is c d v [ln(a(l - (3 2 )/R(v))]] = 



(A41) 
(A42) 



Eq. (|A40[) shows that 3[-F+] depends linearly on s. Consequently, the crack front wave solution s c is a solution of 
Eq. (|A41[) while Eq. (|A42[) is an identity which is satisfied for every real s and v. This result can be checked through 
direct computation of the integral in (|A42j) . Moreover, one has 



^ fell - 3%330)] 



1 S V 1 SI) 



2 (c d a) 2 2 {c s py 
Combining this result with Eq. (|A42|) one gets 

Z{q n (s)} = S -d v [HP(l~p 2 )/R(v))} . 



(A43) 



(A44) 



In the following, we present the computation of G±(s) following the same steps as for the computation of F + [Sc 
Let us introduce the following change of variables in the above defined integrals 



J 



s - vC \ = (s - O 2 



We then have 



and 



leading to 



and 



T(J) = 



7 V[(2 - jjdf - y 1 - j/cj^—fM 

R(v)J(J/c% - 1) 



-sv±iy/ J(J - H 2 ) 
J ~ v 2 



sv i [H 2 {J + v 2 ) - 2v 2 J] 



dJ {J-v 2 ) 2 2^/J(J-H 2 )(J-v 2 ) 2 



± 



(A45) 



(A46) 



(A47) 



(A48) 



(A49) 



Note that we search for solutions for which v 2 < H 2 < c 2 R . Thus, the branch cut runs from c s to Cd and the 
contribution of lnT( J) to the contour integral does not come from the denominator of T(J) which is positive. Then, 
one has 



G±(s) 



1 



tan 



Ay/1 - J/C^J/C 2 ~ 1 



(2 - J /el) 2 



and thus 



G+{s) + G-(s) 



/ dm 


[ £-(■/) " 


tan 1 




/ c 2 





yi - j/cW J/c 2 s - 1 



(2-J/c 2 ) 2 



(A50) 



(A51) 



9 



with 



25R 



[2v 2 (J + s 2 ) - H 2 (J + v 2 )] /(J - v 2 ) [2v 2 J - H 2 (J + v 2 )] /^/J(J-H 2 ) 



(A52) 



£-(J)-i 



{sv) 2 + (^/J{J-H 2 ) + J-v 2 ) 2 (sv) 2 + { V / J(J - H 2 ) + J -v 2 ) 2 



With these algebraic manipulations and simplifications at hand, the numerical computation of 3?[g 33 (s)] and 3?[<j 11 (s)] 
is straightforward. 
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